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ABSTRACT 


During  the  contract  period  the  background  necessary  to  proceed  with  the 

(3*) 

assembly  and  testing  of  the  precision  Hq^l'nuclear  gyroscope  was  completed. 

(3':  4 , 

This  background  work  included  experiments  on  He^He^'  liquid  and  gas  mixtures 

in  a  prototype  apparatus  which  was  modified  to  provide  additional  information 

(3 )  3 ' 

useful  to  the  Het,  gyroscope  research  program.  The  precision  quartz  He^  gyro¬ 
scope  was  designed,  constructed  and  delivered  during  the  period  as  were  the 

components  of  the  airlock  and  cryostat  probe  assembly.  A  new  ultra-low 

-8 

magnetic  field  shield  was  made  which  achieved  2  x  10  G  over  the  volume 

necessary  for  the  Hd_  gyroscope.  This  not  only  exceeds  the  nominal  require- 
3 

ments  of  the  He.,  gyroscope,  but  is  the  lowest  magnetic  field  region  ever  made. 
An  engineering  Ph.D.  thesis  was  completed  by  Captain  Gerald  Shaw,  working  with 
Professor  Daniel  DeBra,  on  a  theoretical  analysis  of  the  cross-axis  reponse 
of  a  three-axis  Hd^  gyroscope  and  on  kinematic  rectification  in  a  nuclear  gyro. 
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I.  Introduction  and  Background 

The  department  of  physics  and  aeronautics-astronautics  of  Stanford 

3 

University  have  undertaken  the  development  of  a  uniquely  stable  He  nuclear 
gyroscope  based  on  low-temperature  technology.  This  development  effort 

makes  use  of  experience  we  have  already  obtained  in  techniques  of  polarizing 

3  3 

He  by  optical  pumping,  obtaining  long  He  nuclear  relaxation  times,  obtain¬ 
ing  low-magnetic-field  regions  in  superconducting  shields  and  selecting 
magnetically  clean  materials  for  use  in  these  low-field  environments,  test¬ 
ing  and  working  with  precision  fused  quartz  gyro  assemblies  and  using  SQUID 
magnetometry  for  gyro  readout.  Funding  for  the  previous  work  has  been 
provided  on  large  part  by  AFOSR  and  NASA. 

Prior  to  the  current  AFOSR  funding,  a  crude  prototype  of  the  cryogenic 

3  3  3  4 

He  nuclear  gyroscope  (He  NG)  was  built.  Using  this  apparatus.  He  -He  mix¬ 
tures  were  condensed  into  a  4°  K  low-field  region  and  the  resulting  processing 
magnetization  was  observed  with  a  signal-to-noise  ratio  in  excess  of  103  by 
use  of  a  SQUID  magnetometer.  Although  this  prototype  was  too  flawed  (in  large 

part  due  to  magnetic  contaminants)  to  function  as  a  gyro,  it  was  found  that 
3  4 

0.07%  He  -  liquid  He  mixture  had  a  relaxation  time  of  40  hours  in  a  1  cm 

diameter  Pyrex  cell  when  the  effect  of  magnetic-gradient-induced  relaxation 

was  deducted.  This  result  was  dominated  by  wall-induced  relaxation  since 

the  subsequent  use  of  a  solid  H^  wall  coating  increased  this  relaxation  time 

to  more  than  five  days  (140  hours).  This  last  result  is  encouragingly  close 

to  the  theoretically  estimated  maximum  relaxation  time  under  these  conditions 

3  4 

of  170  hours  due  to  He  -He  dipolar  interactions. 
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In  addition  to  these  experimental  measurements,  a  general  analysis  of 
the  characteristics  of  a  He~  NG  based  on  current  low-temperature  technology 
was  made  for  the  purpose  of  evaluating  its  potential  for  use  in  a  future 
nuclear  electric  dipole  moment  experiment.  This  analysis  considered  two 
possible  modes  of  operation  of  the  He  NG  and  the  results  of  one  of  these 
analyses  are  presented  in  Section  III.  The  details  of  both  the  theoretical 
analyses  and  the  relaxation  time  measurements  can  be  found  in  the  Ph.D. 
thesis  of  M.  A.  Taber.  A  summary  of  the  experimental  results  was  also 
presented  at  the  15th  International  Conference  on  Low  Temperature  Physics 
and  published  in  the  Proceedings. 

In  1978,  work  was  begun  with  support  from  AFOSR  on  the  development  of 

4 

a  working  precision  cryogenic  He  NG.  This  work  has  resulted  in  a  gyro  housing 
design  (Fig.  1).  This  design  and  its  tolerance  sp  wifications  are  the  result 
of  the  elaboration  of  our  previous  general  analysis  plus  additional  analytical 
work.  Because  of  our  desire  to  investigate  the  possibility  of  ultimately 
using  this  instrument  for  a  demanding  physics  experiment,  design  emphasis 
has  been  placed  on  a)  magnetic  field  stability,  b)  long  relaxation  time, 
c)  good  readout  resolution,  and  d)  flexibility  anf  ease  of  modification. 

The  housing  was  designed  in  collaboration  with  Mr.  Don  Davidson  of  Optical 
Instrument  Design,  Inc.  who  supervised  the  fabrication  of  the  unit. 
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Fig.  1.  Schematic  cross  sectional  view  of  the  He  gyro  housing.  Coil 
forms  Z^  -  Z^  are  axially  symmetric  about  the  Z  axis  and  are  configured 
to  generate  the  primary  6th-order  field.  Coils  X^,  X2(not  shown)  and 
Y2  are  symmetric  about  the  X  and  Y  axes,  respectively,  and  provide 
4th-order  fields  for  correction  and  initialization  purposes.  Coil  P  is 
the  readout  coil.  Not  9hown  is  the  valve  assembly  which  resides  in  the 
fill-line  adapter.  All  parts  shown  are  fused  silica. 


II.  Description  of  the  Low  Temperature  He3  Nuclear  Gyro 

It  is  a  well-known  result  of  quantum  mechanics  that  the  equations 
of  motion  of  the  expectation  value  of  a  nuclear  magnetic  dipole  moment 
U  =  Y»I,  where  nl  is  the  nuclear  spin  vector  operator,  are 

d<p>/dt  =  y<M>  x  B  [1] 

1 

when  the  nucleus  is  placed  in  a  magnetic  field  B.  This  of  course,  is 
formally  identical  to  the  classical  equations  of  motion  that  would  be 
obeyed  by  any  freely  spinning  body  that  possessed  a  magnetic  dipole 
moment  aligned  along  its  angular  momentum  vector.  Another  basic  theor¬ 
etical  result  is  that  in  the  case  of  spin-1/2  particles  (such  as  the  He3 

2 

nucleus),  there  can  be  no  higher  multipole  moment  than  the  dipole. “  Thus 
[1]  represents  the  only  external  interaction  that  is  necessary  for  us  to 
consider. 

3 

If  a  large  number  of  non-interacting  He  atoms  are  collected  in  a 
sample  volume  and  are  subjected  to  a  homogeneous  magnetic  field,  then 
provided  the  components  of  the  total  angular  momentum  are  large  compared 
to  h,  the  motion  of  the  total  angular  momentum  can  be  treated  classically. 
Thus  we  may  replace  the  expectation  values  of  the  individual  spin  vectors 
in  [1]  with  the  classical  variable  M,  the  sample  magnetization: 

dM/dt  =  yCM  x  B).  [2] 

Equation  [2]  raises  an  intriguing  question:  is  it  possible  to  use 

3 

such  a  collection  of  He  atoms  to  make  a  practical  gyroscope?  In  order 


3  3 

to  make  a  useful  He  nuclear  gyroscope  (He  NG),  we  must  at  least  approxi¬ 
mate  the  various  conditions  that  were  noted  in  arriving  at  equation  [2], 
and  additionally  the  magnetic  environment  must  be  strictly  controlled.  In 

3 

practical  terms,  this  means  the  following  requirements  must  be  met:  (1) 

A  sample  containing  a  large  number  of  polarized,  rapidly-diffusing  He3 
atoms  must  be  obtained.  (2)  The  components  of  the  magnetization  must  have 
usefully  long  relaxation  times.  (3)  The  polarized  He3  must  be  maintained 
in  a  spherical  sample  cell,  and  any  nearby  materials  that  have  large  mag¬ 
netic  susceptibilities  (paramagnetic  or  diamagnetic)  must  also  have  at 
least  cubic  symmetry.  (4)  If  there  is  to  be  relative  motion  between  the 
sample  magnetization  and  the  gyro  case,  materials  that  can  cause  signifi¬ 
cant  radiation  damping  must  not  be  in  close  proximity  to  the  sample.  (5) 

A  sensitive  means  of  measuring  the  direction  of  the  sample  magnetization 
without  perturbing  the  magnetization  to  an  unacceptable  degree  must  be 
provided.  And  finally,  (6)  the  average  magnetic  field  over  the  sample 
volume  (exclusive  of  that  due  to  the  sample  itself)  must  either  be  intrin¬ 
sically  stable  or  independently  measured  and  controlled  with  sufficient 
stability.  In  view  of  the  relaxation  time  requirement  in  (2),  it  is  also 
important  that  the  magnetic  field  gradients  be  kept  small. 

3 

There  are  at  least  two  ways  in  which  a  He  NG  might  be  formulated. 

One  way  would  be  to  reduce  the  average  magnetic  field  over  the  sample 
volume  to  zero.  In  this  case,  it  can  be  seen  from  equation  [2]  that  M 
should  remain  fixed  in  direction,  and  rotation  of  the  gyro  housing  about 
any  axis  orthogonal  to  M  could  be  measured.  Although  this  approach  has 
the  advantage  of  two-axis  sensitivity,  an  analysis  of  its  potential  per¬ 
formance  (as  implemented  by  low  temperature  techniques)  indicates  that  in 
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4 

most  respects  it  will  be  inferior  to  the  alternative  approach. 

The  alternative  approach  is  to  apply  a  very  constant,  uniform  mag¬ 
netic  field  B  =  B  R.  By  letting  u>  =  -yB  (y  =  -2.04  x  10*  rad  sec'1  G  1 
o  o  o  o 

for  He^)  equation  [2]  can  be  written  in  this  case 

dM/dt  =  -  (M  *  uiQ)  .  [2 '  J 

Equation  [2']  is  readily  solved  by  the  classical  technique  of  taking  the 
time  derivative  of  M  in  another  frame  that  is  rotating  with  respect  to 
the  inertial  frame  with  angular  velocity  ft: 

(^“'’inertial  '  ^’rot  *  5  *  3  • 

thus 

(dS/dt)rot  =  M  x  (ft  -  Zq) .  [3] 

Equation  [3]  shows  that  M  is  stationary  in  the  frame  that  is  rotating  with 

angular  velocity  u)q  irrespective  of  the  orientation  of  M  with  respect  to 

B  .  (Of  course,  M  is  a  constant  of  the  motion.)  Thus  M  precesses  about 
o  z 

the  z  axis  at  the  Larraor  precession  frequency  u)q  with  respect  to  inertial 
space . 

We  now  consider  the  field  sources  of  BQto  be  fixed  with  respect  to 
the  gyro  case  and  ask  what  happens  when  the  gyro  case  is  rotated  with  re¬ 
spect  to  the  inertial  frame.  If  the  case  is  rotated  at  a  rate  ft  =  ftK  we 
see  by  comparing  £ 2 * ]  and  [3]  that  the  precession  frequency,  as  measured 
in  the  frame  of  the  gyro  case  is  shifted  by  an  amount 

<$u)  =  ft  .  [4] 

o  1  J 
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By  integrating  [4]  with  respect  to  time  we  also  see  that  if  the  case  is 

rotated  about  the  z  axis  by  an  angle  A<f>,  then  the  Larmor  phase  will  also 

shift  the  same  amount.  Thus,  by  initially  establishing  M  in  the  x-y 

plane  and  continually  monitoring  the  phase  of  M  (or  M  ) ,  rotations 

x  y 

about  the  z  axis  can  be  measured.  It  is  clear,  however,  that  since  wo 
is  not  fixed  by  any  constant  of  nature,  its  phase  and  frequency  must  be 
initially  established  by  measurement  when  the  gyro  case  is  mounted  on  an 
inertial  platform  or  is  at  least  subjected  to  an  accurately  known  rate  of 
rotation. 

Now  if  the  case  is  rotated  about  the  x  or  y  axis  with  a  velocity 
such  that  |^|  «  |uj  |,  then  it  is  seen  from  [3]  that  in  the  frame  of  the 
gyro  case,  M  will  precess  about  toQ  -  d  at  a  rate  wq[1  +  0(£2/<i>o)^]  •  Thus 
M  will  "adiabatically  follow"1  Bq  with  only  a  second  order  shift  in  the 
Larmor  frequency.  This  approach  therefore  yields  a  single-axis  gyroscope 
which  is  limited  to  rotation  rates  that  are  small  compared  to  u)Q.  It  is 
this  type  of  gyroscope  that  we  will  now  discuss  in  greater  detail  within 
the  context  of  low  temperature  techniques. 

The  low  temperature  He^HG  is  depicted  schematically  in  Fig.  2.  Al¬ 
though  low  temperature  techniques  introduce  certain  complications,  they 
offer  two  significant  advantages  to  the  He^NG.  The  first  is  that  the 
flux-quantization  property  of  type  I  superconductors  means  that  the  applied 
field  is  intrinsically  stable  except  for  two  effects  that  will  be  discus¬ 
sed  below:  the  variation  of  dimensions  and  penetration  depth  with  tempera¬ 
ture.  This  stability  makes  it  unnecessary  to  independently  monitor  the 
magnetic  field. 

The  other  advantage  is  that  a  SQUID  (Superconducting  Quantum  Inter - 
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SUPERCONDUCTING 
COUPLING  CIRCUIT 


Fig.  2. 


Schematic  of 


He  nuclear  gyroscope 


ference  Device)  magnetometer  can  be  used  for  readout.  This  device  can 
be  coupled  to  the  sample  magnetization  with  a  superconducting  circuit. 

The  output  of  the  magnetometer  is  fed  back  into  this  circuit  so  that 
the  reaction  currents  that  would  normally  flow  due  to  motion  of  the 
sample  magnetization  are  reduced  by  a  factor  equal  to  the  loop  gain  of 
the  magnetometer.  This  readout  technique,  then,  utilizes  the  high  sensi¬ 
tivity  of  the  SQUID  for  accurate  readout  and  at  the  same  time  minimizes 
the  effect  of  the  coupling  circuit  on  the  sample  magnetization. 

For  the  purposes  of  discussion  we  will  now  describe  a  possible  ver- 

3 

sion  of  the  low  temperature  He  NG  in  some  detail.  The  values  of  the 
various  parameters  that  are  specified  here  form  the  bases  for  the  esti¬ 
mates  of  the  gyro  performance  that  will  be  quoted  later.  It  should  be 
noted  that  the  choices  of  parameter  values  were  based  on  (a)  reasonable 
compromises  between  various  tradeoffs,  although  certain  applications  may 
require  different  compromises  than  those  selected  here;  (b)  capabilities 
of  our  current  optical  pumping  and  low  temperature  laboratory  equipment; 
and  (c)  laboratory  experience  using  techniques  that  are  similar  (if  not 
identical)  to  those  that  are  proposed  here. 

The  He^NG  that  is  being  contemplated  consists  of  a  fused  quartz 
housing  with  a  3.8  cm  diameter  sample  cell  that  has  been  tumble  lapped 
and  assembled  to  maintain  a  sphericity  of  +  0.13  pm.  Permanent  leak- 
tight  assembly  is  achieved  by  optical  contacting.  The  housing  includes 
as  integral  elements  the  pickup  loop  for  the  magnetometer,  field  coils, 
and  a  15  cm  diameter  spherical  superconducting  shield  on  the  outer  sur¬ 
face.  Because  all  of  the  conductors  in  the  quartz  housing  are  to  be 
superconducting  and  therefore  perfectly  diamagnetic,  they  must  either  be 
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physically  small  or  arranged  to  maintain  cubic  symmetry  if  excessive  pre¬ 
cession  of  the  sample  magnetization  is  to  be  avoided. 

The  quartz  housing  is  to  be  located  in  a  cryostat  probe  that  can  be 
cooled  to  4.2  K  by  lowering  into  a  dewar  that  is  lined  with  a  degaussed 
superconducting  lead-foil  shield.  This  shield  is  degaussed  by  use  of  an 
iterative  expansion  technique  that  has  been  perfected  in  our  laboratory,5 

and  can  be  expected  to  provide  a  low-field  region  with  a  magnetic  field 
-8 

of  -  5  x  10  G.  It  is  estimated  that  this  field  should  have  a  gradient 
-8-1 

on  the  order  of  10  G  cm  in  the  center.  If  the  spherical  gyro  shield 
is  carefully  cooled  in  this  ambient  field,  it  will  trap  a  field  of  similar 
magnitude  and  configuration. 

In  our  analysis,  we  have  assumed  a  value  for  the  Larmor  frequency  of 
o)o  =  1  rad  sec'1  (Bq  =  5  *  10'5  G) .  This  number  is  primarily  constrained 
by  magnetic  gradient -induced  spin  relaxation  considerations  (see  below) . 

It  would  be  possible  to  increase  u)q  significantly  with  no  degradation  of 
gyro  performance  or  reduction  of  dimensional  tolerance  requirements  if 
the  diameter  of  the  gyro  housing  (i.e.,  the  spherical  shield)  were  made 
larger. 

The  He5  can  be  polarized  at  room  temperature  using  optical  pumping*3 

and  admitted  to  the  sample  cell  via  a  glass  filling  tube.  A  mechanical 

closure  in  the  form  of  a  tapered  quartz  plug  can  be  opened  to  admit  the 

sample  and  then  closed  to  prevent  diffusion  in  and  out  of  the  sample  cell. 

3  4 

Based  on  our  experience  with  the  condensation  of  polarized  He  -  He 
mixtures,  it  is  estimated  that  if  the  He5  is  polarized  to  10%  at  room 
temperature  it  will  have  -  2%  initial  polarization  in  the  sample  cell. 

If  it  is  required  that  the  initial  internal  magnetic-flux  density  be 

3  x  10  G,  then  a  He-3  density  of  3  x  10  mol  cm  is  necessary  in  the  sample 

cell. 

10 


It  will  undoubtedly  be  essential  that  a  cryogenic  wall  coating  of 

solid  be  used  in  order  to  reduce  wall-induced  nuclear  relaxation  of 

37  4 

the  He  .  It  is  possible  that  a  certain  amount  of  He  gas  could  be  used 

to  reduce  further  wall-induced  relaxation  as  will  be  discussed  below. 

The  last  element  that  must  be  specified  in  order  to  estimate  the 
3 

He  NG  performance  are  the  noise  and  drift  characteristics  of  the  SQUID 
For  this  purpose  we  have  used  the  Clarke,  Goubau,  and  Ketchen  (CGK)  data 

g 

on  dc  tunnel- junction  SQUIDs.  It  should  be  noted,  however,  that  the  dc 
tunnel-junction  SQUIDs  are  still  only  laboratory  devices  and  are  not  yet 
commerically  available.  Currently  we  are  using  rf-biased  hybrid  SQUIDs 
which  have  a  noise  level  approximately  three  times  higher  than  the  CGK 
result  when  referred  to  the  input. 

*  *  T 

Having  described  a  possible  low  temperature  He  NG,  we  will  now  sum¬ 
marize  the  results  of  our  estimates  of  three  important  performance  para- 
meters  based  on  this  description:  (1)  He  nuclear  relaxation  time,  (2) 
readout  resolution,  and  (3)  stability  of  the  Larmor  frequency,  i.e.,  drift 
and  precession  effects. 

1.  He^  Nuclear  Relaxation  Time 

The  relevant  relaxation  time  is  T^»  the  transverse  relaxation  time. 
There  are  three  significant  relaxation  mechanisms:  intrinsic  relaxation, 
wall-induced  relaxation,  and  relaxation  due  to  magnetic  field  gradients. 

In  the  first  two  cases  T^  =  T  ,  the  longitudinal  relaxation  time.  In  the 
case  of  relaxation  due  to  field  gradients,  however,  T^  and  are  not 
generally  equal,  although  in  the  present  case  they  should  be  of  the  same 
order  of  magnitude. 

Both  intrinsic  and  gradient -induced  relaxation  are  relatively  well- 

understood.  In  the  case  of  intrinsic  relaxation,  the  relaxation  rate  of 

„  3  9 

pure  He  gas  is  given  approximately  by 
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[5] 


T"1  =  2y^2ICI  ♦  l)(7r&n)*n3<f2 

-1  3 

where  I  =  %  is  the  nuclear  spin  quantum  number,  0  =  (kT)  ,  n^  is  the  He 
number  density,  and  d  is  the  atomic  hard-core  diameter.  We  note  that  of 
the  controllable  parameters  in  equation  [5],  0  and  n^,  is  the  most 
important.  Under  the  proposed  conditions  for  the  He^NG,  [5]  yields  the 
estimate  Tj  =  7  x  io6  sec  =  80  d. 

The  expressions  for  nuclear  relaxation  due  to  diffusion  through  mag¬ 
netic  field  gradients^  are  not  particularly  difficult  to  obtain  but  are 
quite  cumbersome  for  an  arbitrary  field  configuration.  For  simplicity, 
we  will  limit  the  present  discussion  to  the  case  of  a  uniform  field  plus 
an  axially  symmetric  uniform  gradient:  B(r)  =  BqK  +  g(2ix  -  jy  -  Jcz)  . 

In  this  case 

T'J  <  O.ISyYrV1  [6J 

3  -1 

where  D  is  the  He  diffusion  coefficient  and  T.„  is  the  transverse  relaxa- 

3  2 

tion  rate  due  to  the  gradient  g.  For  the  proposed  He  density  D  :  1  cm 

sec  *,  and  we  estimate  *  IO2  sec.  This  is  obtained  with  the  assump- 
-8  - 1 

tion  that  g  *  10  G  cm  .  It  is  important  that  care  be  taken  to  avoid 

114 

ferromagnetic  contaminants  near  the  sample  volume. 

In  the  case  of  wall-induced  relaxation,  the  situation  is  not  so 

certain.  The  only  data  that  appear  to  be  available  of  He"*  relaxation 

measurements  with  a  solid  Hj  wall  coating  are  the  measurements  of 

Barbe,  Laloe,  and  Brossel7  (BLB) ,  and  our  Tj  measurements.4  The  BLB 

3  -2 

measurements  were  made  with  a  He  density  ~  10  *  that  being  proposed 

here.  They  obtained  a  relaxation  time  (T^)  of  60  hours  in  a  3  cm  diameter 
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cell  that  was  completely  dominated  by  wall-induced  relaxation.  Our  results 

on  the  other  hand,  were  obtained  with  a  He^  density  ~  7*  that  being  present 

4 

ly  considered  and  in  the  presence  of  liquid  He  as  a  buffer.  We  obtained 

a  relaxation  time  of  140  hours  (in  a  1  cm  diameter  cell)  that  compares  with 

a  (crude)  theoretically-estimated  intrinsic  relaxation  time  of  ~  170  hours. 

12 

From  these  and  other  data,  we  conclude  (a)  the  data  are  presently 

l 

insufficient  to  predict  with  accuracy  what  the  wall-induced  relaxation 

rate  (and  thus  the  resultant  relaxation  time)  will  be  under  the  conditions 
3  4 

proposed  for  the  He  NG,  and  (b)  a  He  gas  buffer  could  be  used  to  reduce 
further  wall-induced  relaxation  if  it  should  be  necessary.  Unfortunately, 

4 

unlike  the  H^  wall  coating.  He  has  a  liability  as  well  as  an  advantage: 

4  3 

the  addition  of  He  decreases  the  He  diffusion  coefficient,  and  the 
gradient - induced  relaxation  rate  is  increased  proportionately.10  Thus, 

4 

optimal  He  density  depends  on  the  relative  potency  of  the  gradient-  and 

wall-induced  relaxation  mechanisms.  We  estimate,  however,  that  with  a 
-8-1  4 

gradient  of  10  G  cm  ,  the  optimal  He  density  will  lead  to  a  relaxation 
time  in  excess  of  105  sec  but  probably  not  more  than  106  sec. 

2.  Readout  Resolution 

Since  the  desired  information  is  encoded  at  a  non-zero  frequency, 

v  =o)  /2ir,  we  can  estimate  the  readout  resolution  in  the  limit  of  a 
o  o 

large  signal-to-noise  ratio  with 


%  *tVVMt>  ^  “"o’'  I?1 

where  At  is  the  measurement  time,  B  is  magnitude  of  the  internal  mag- 

mo 

netic  flux  density:  B  =  (8u/3)M,  and  SD(v  )  is  the  one-sided  magneto- 

mo  D  O 

meter  noise  spectral  density  referred  to  BmQ.  This  expression  is  based  on 
the  assumption  that  VQAt  >>  1.  The  necessity  for  this  assumption  derives 


from  the  fact  that  this  device  is  blind  to  small  rotations  when  the  instan¬ 
taneous  orientation  of  the  magnetization  is  orthogonal  to  the  pickup  coil. 
This  assumption  is  unnecessary  and  the  expression  in  (  7]  is  reduced  by  a 
factor  of  /2  if  two  orthogonal  pickup  coils  and  SQUIDs  are  utilized  instead 
of  one.  If  this  were  the  case  (as  it  undoubtedly  would  be  for  a  practical 
instrument),  the  only  bandwidth  limitation  would  be  that  imposed  by  the 
SQUID  electronics. 

If  we  assume  that  the  pickup  coil  consists  of  a  compact  filamentary 

L  ^ 

loop  in  the  y-z  plane,  then  is  given  by 

S*  *  (2/it)  (2En<)  V/2r’3  , 

where  E^fv)  is  the  magnetometer  energy  sensitivity  at  frequency  v,  a  is 

the  pickup  coil  radius,  Rq  is  the  sample  cell  radius,  and  <  is  a  slowly 

varying  factor  relating  the  coil  cross  section  to  its  inductance.  The  CGK 

data  show  that  Er  is  approximately  white  and  equal  to  h  10-30JHz  1  at  vq 

and  we  estimate  *  3  *  10  11  G  Hz~^.  Thus  a,  -  10  5 (At)  ,  if  B  »  3 
b  <p  mo 

x  10  ^  G.  As  the  magnetization  decays  with  the  relaxation  time  T^,  the 
readout  resolution  will  also  deteriorate. 

3.  Stability  of  the  Larmor  Frequency 

The  most  important  sources  of  variation  in  the  Larmor  frequency  are 
likely  to  be  (a)  variations  in  Bq  due  to  mechanical  motion  arising  from 
external  perturbations  or  dimensional  or  penetration-depth  variation  arising 
from  thermal  fluctuations,  (b)  the  varying  magnetic  field  generated  in 
the  pickup  loop  because  of  noise  currents  injected  into  the  input  circuit 
by  the  magnetometer,  (c)  interaction  of  the  sample  magnetization  with  it¬ 
self  because  of  asymmetries  in  the  sample  cell,  superconducting  parts 
(including  the  shield)  and  because  of  finite  magnetometer  loop  gain,  and 
(d)  the  London  moment.  It  is  easiest  to  consider  these  effects  separately. 
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a.  Effect  of  Dimensional  and  Temperature  Stability  on 

Although  a  detailed  design  of  the  gyro  case,  field  coils,  etc.  is 

needed  to  estimate  the  magnitude  of  dimensional  variations  and  their 

effect  on  the  Larmor  frequency,  we  can  make  some  brief  remarks.  First, 

because  the  magnetic  field  gradient  must  be  small  because  of  relaxation 

considerations,  the  effect  of  small  displacements  of  the  sample  cell  with 

respect  to  the  spherical  shield  will  be  minor.  For  instance,  with  a 
-8  - 1 

gradient  of  10  G  cm  ,  a  relative  displacement  as  large  as  0.1  pm  would 

-9  -1 

produce  a  frequency  shift  of  only  6u)q  -  2  x  10  rad  sec  .  Secondly, 

-8 

the  linear  expansion  coefficient  of  fused  silica  at  4.2  K  is  a  :  -3  *  10 
K  Now  if  the  field  coils  were  to  consist  of  a  persistent  supercon¬ 

ducting  circuit  that  is  deposited  on  a  quartz  surface,  then  6Bq/Bo  = 

6u>o/o)  =  -2a<$T  =  6  x  10  8(6T),  where  6T  is  in  degrees  Kelvin. 

There  is  an  additional  effect  that  can  cause  B  to  vary  as  a  function 

o 

of  temperature:  the  penetration  depth  (i.e.,  the  average  depth  that  a 
magnetic  field  penetrates  the  surface  of  a  superconductor)  has  a  weak  tem¬ 
perature  dependence  for  temperatures  T  «  T  ,  Tc  being  the  transition  tem¬ 
perature  of  the  superconductor.  This  situation  is  most  simply  modeled 
by  considering  the  average  magnetic  field  produced  by  a  persistent  current 
in  a  superconducting  ring  that  is  thick  compared  to  the  penetration  depth. 
The  model  then  takes  the  effective  inner  radius  to  be  equal  to  the  inner 
radius  of  the  ring  plus  the  penetration  depth.  Since  the  flux  through  the 
ring  is  conserved  we  have 


26r  2  fdX-\, 

~ =  ■  r  qrj6T  ■ 


where  r  is  the  inner  radius  of  the  ring,  and  A  =A(T)  is  the  temperature- 
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dependent  penetration  depth,  which  for  our  purposes  is  given  with  sufficient 

.  I4 

accuracy  by 

A(T)  =  A0[l  -  (T/Tc)4]"*  . 

In  the  limit  that  T  «  T  , 


If  the  shield  is  made  out  of  niobium,  Tc  =  9.46  K,  and  XQ  =  4.7  x  10”6  cm.15 
Taking  r  =  6  cm  we  then  estimate 

6B  /B  -  -  3  x  10'8(6T) 

OO 

for  T  =  4.2  K.  This  is  of  the  opposite  sign  and  approximately  one-half 

the  size  of  the  effect  due  to  the  coefficient  of  expansion  of  the  quartz. 

Temperature  regulation  at  the  level  of  +  10  mK  will  therefore  lead  to 

variations  in  o>  on  the  order  of  3  x  10-I°rad  sec-1, 
o 

b.  Effect  of  Magnetometer  Noise  and  Drift  on  o>o 

As  we  have  noted  previously,  it  is  advantageous  to  include  the  super¬ 
conducting  coupling  circuit  as  an  element  in  the  magnetometer  feedback  loop 
in  order  to  reduce  reaction  currents.  In  doing  this,  the  noise  currents 
injected  in  the  coupling  circuit  by  magnetometer  noise  and  drift  become 
accurately  known  (in  the  limit  of  large  loop  gain)  and  we  can  readily 
estimate  their  effect  from  published  magnetometer  data. 

Since  the  normal  to  the  plane  of  the  pickup  loop  should  be  orthogonal 
to  Bq,  the  magnetic  field  generated  by  small  currents  in  the  pickup  loop 
should  ideally  have  only  a  second  order  effect  on  |Bq| .  Because  of  con¬ 
struction  tolerances  and  because  of  the  ambient  field  trapped  in  the  shield, 
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however,  the  unit  normal  of  the  pickup  loop  will  have  a  small  component 

in  the  direction  of  B  of  the  order  of  10  .  This  will  be  much  more  im- 

o 

portant  than  the  second  order  effect. 

An  analysis  based  on  the  CGK  data  indicates  that  the  effect  of  linear 
magnetometer  drift  becomes  dominant  over  that  of  random  noise  after  about 
an  hour  or  so.  If  we  now  use  the  average  drift  rate  over  a  20  hour  period 
reported  by  CGK,^  we  estimate  that  the  Larmor  frequency  will  have  changed 
by  '  7  x  io  ^  rad  sec  *  after  20  hours.  This  means  that  the  average  gyro 
drift  rate  over  a  20  hour  period  will  be  ~4  *  10  rad  sec  1 . 

It  should  be  pointed  out  that  the  true  gyro  error  is  significantly 
less  than  this.  Because  the  signal  phase  is  the  integral  of  the  Larmor 
frequency,  only  low  frequency  noise  components  and  drift  introduce  signi¬ 
ficant  error.  Since  the  magnetometer  loop  gain  is  quite  large  at  low 

5 

frequencies  10  ) ,  it  should  be  possible  to  determine  accurately  the 
relationship  between  u)Q  and  the  low  frequency  noise  and  drift  that  is  seen 
at  the  magnetometer  output.  Once  this  is  done,  errors  introduced  by  this 
source  can  be  reduced  even  further. 

c.  Effect  of  Gyro  Housing  Asymmetries  on  u)q 
It  can  he  shown  that  the  effect  of  asymmetries  in  the  sample  cell 
and  in  superconducting  components  causes  a  frequency  shift  that  is  pro- 
portional  to  yMz  =  y|M(cos0.  The  fact  that  the  magnetometer  loop  gain 
is  finite  has  a  similar  effect.  This  leads  to  a  slow  variation  in  w 

o 

since  Mz  decays  with  a  time  constant  T^.  As  an  example  of  this  type  of 
effect  we  will  write  down  the  appropriate  expression  for  the  effect  of 
sample  cell  asymmetry.  It  is  estimated  that  other  contributions  will 
be  of  the  same  order  or  smaller. 
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If  the  radius  of  the  sample  cell  is  expressed  in  the  form 


R(M)  =  Rq[1 


I/O 


where  «  1,  8^  ant*  are  spherical  harmonics  (the 

polar  axis  being  in  the  direction  of  Bq) ,  then  it  can  be  shown  that  only 
8„rt  has  an  effect  in  first  order.  By  defining 

4.U 

{R20  ’  Ro620Y20<6=0>  ’  '/573?  R„S20 

4 

the  frequency  shift  due  to  sample  cell  asymmetry  is  simply 


«“o.  =  -  (9/10)(«*M/R))YB.ocoSe  '  I*) 

In  order  to  evaluate  the  gyro  drift  from  this  source  it  is  necessary 
to  know  cos9  and  its  evolution  in  time.  It  is  reasonable  to  assume  that 
|cos0Q]  <  10  when  the  gyro  is  initialized.  There  are,  however,  three 
mechanisms  which  can  cause  the  magnetization  to  move  away  from  the  x-y 
plane:  (a)  radiation  damping  due  to  the  magnetonjeter  gain  being  finite 
and  having  a  non-zero  phase  shift,  (b)  the  component  of  the  magnetic-noise 
field  at  frequency  u>o/2tt  that  originates  from  the  magnetometer  feedback, 
and  (c)  the  component  of  random  angular  displacements  of  the  gyro  case 
about  the  x  and  y  axes  at  frequency  ujo/ 2tt  .  The  first  of  these  mechanisms 
will  cause  (0  -  tt/2)  to  vary  linearly  in  time,  whereas  the  latter  two 
will  cause  0  to  execute  a  random  walk  so  that  the  variance  of  (0  -  tt/2) 
will  increase  linearly  in  time. 

Under  the  conditions  that  the  magnetometer  loop  gain  magnitude  is 
5  -2 

~  10  with  an  argument  £  10  rad  at  o>o/ 2tt ,  and  the  platform  angular 
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-6  -h 

noise  spectral  density  is  £  10  rad  Hz  7  ( also  at  uq/2it),  then  it  is 

estimated  that  0  will  not  change  more  than  10~3  rad  over  a  period  of 

10^  sec  (12  d) .  Since  these  conditions  are  quite  reasonable,  we  will 

assume  that  |cos0|  <  2  x  10  in  equation  [8].  By  taking  6R7q  *  0.13  pm, 

R  ■  1.9  cm,  B  =  3  x  io~**  G,  we  then  obtain  the  estimate  |6o)  I  <  10"® 
o  mo  1  oa1 


rad  sec 


-1 


d.  Effect  of  the  London  Moment 

3 

If  a  spherical  superconductor  such  as  the  shield  for  the  He  NG  is 

rotated  with  an  angular  velocity  then  a  uniform  field  will  appear  in 

18 

the  interior  according  to  the  expression  derived  by  London: 


Here,  m  and  e  refer  to  the  mass  and  charge  of  the  electron,  and  the  numeri 
cal  value  of  the  coefficient  is  1.14  x  10~7  G  rad'1  sec.  Writing  = 

-  YBl,  we  have  the  frequency  shift  due  to  the  London  moment: 

y(—  )$  =  2.6  x  10'3  S  . 

ol  v  e  ' 


Even  though  this  effect  is  large,  it  does  not  represent  a  true  gyro  error 
since  it  is  a  well-known  systematic  effect  that  can  be  readily  corrected 
in  the  data  processing. 

3 

In  order  to  make  use  of  the  estimated  stability  of  the  He  NG,  it  is 
necessary  that  the  stability  of  the  clock  used  to  measure  the  phase  of 
the  He3NG  output  signal  not  be  a  limiting  factor.  There  are,  however,  a 

number  of  clocks  that  have  fractional  frequency  stability  (Allan  variance) 

-11  19 

of  10  or  better  for  sampling  times  in  excess  of  1  sec.  Thus  there 


should  not  be  any  problem  in  principle,  but  care  will  be  needed  in  making 


the  phase  measurement  if  the  intrinsic  stability  of  both  the  clock  and  the 
He"*  NG  are  to  be  successfully  utilized. 
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III.  Work  during  Contract  Period 

1.  New  Low  Field  Facility 

3 

One  of  the  most  important  aspects  of  the  cryogenic  He  NG  is  the  ultra-low 
field  facility.  Over  the  past  number  of  years  we  have  developed  considerable 
experience  in  techniques  for  cooling  superconducting  shields  such  that  there 
is  very  little  flux  trapped  in  the  shield.  Recently,  however,  we  have  devised 
an  improved  technique  using  a  new  cooling  apparatus.  During  this  past  period, 
this  technique  has  been  experimentally  tested.  This  new  process  is  designed 
to  significantly  reduce  thermal  gradients  which  can  produce  excess  trapped 
field  during  the  cooling  process.  As  a  result  of  this  development,  we  have 
successfully  produced  an  8"  diameter  shield  for  the  He3  nuclear  gyro  which  has 

_  g 

a  field  level  of  2x10  gauss  over  a  substantial  volume  (Fig.  3).  This  sur- 

_  g 

passes  our  target  of  5x10  gauss  and  to  the  best  of  our  knowledge  is  the 
lowest  absolute  magnetic  fie’d  region  ever  made. 

2.  Precision  Gyro  Housing  and  Cryostat  Probe 

During  this  period  the  precision  ground  fused  silica  components  of  the  new 
gyro  housing  were  received  at  Stanford  (Fig.  4).  The  sample  cell  can  now  be 
tumble  lapped  to  within  5  microinches  of  sphericity  and  assembled  by  optical 
contacting.  The  only  additional  work  that  might  precede  final  assembly  of  the 
housing  components  is  the  deposition  of  the  superconducting  shield  material  on 
the  inner  hemispherical  surfaces  of  the  outer  shell.  This  step,  however,  can 
(and  probably  will)  be  delayed  until  after  initial  tests  have  been  made  on  the 
precision  gyro  assembly. 

In  addition  to  the  gyro  housing,  a  cryostat  probe  that  will  be  used  for 
testing  pruposes  is  under  construction  along  with  an  airlock  assembly.  The 
airlock  assembly  allows  the  cryostat  probe  to  be  slowly  lowered  into  the  dewar 
containing  the  low  field  so  as  to  minimize  thermal  shock  of  the  gyro  components 
and  minimize  boiloff  of  liquid  helium. 


Fig.  4.  Precision-ground  fused-silica  components  to  be  assembled  for  the 
new  gyro  housing. 
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3.  Engineering  Analysis  of  the  Cryogenic  He^NG 

During  this  period  a  final  copy  of  the  engineering  Ph.D.  thesis  of 
Capt.  G.  Shaw  was  submitted  to  the  department  of  aeronautics  and  astronautics. 
This  work  investigated  kinematic  rectification  in  nuclear  gyros  (the  response 
to  soall-angle  coning  of  the  gyro  case  about  the  input  axis),  the  mechanism 
of  cross-coupling  (i.e.,  how  to  reduce  the  error  introduced  by  the  rotation 
of  the  case  about  an  axis  orthogonal  to  the  input  axis) ,  and  gyro  drifts 
induced  by  acceleration  or  thermal  gradients  in  the  presence  of  a  magnetic 
field  gradient.  A  practical  compensation  technique  for  the  coning  effect 
has  been  devised.  These  effects  can  occur  in  any  kind  of  nuclear  gyro. 

Copies  of  the  thesis  are  being  included  as  a  part  of  this  report.  In  addition 
to  the  analysis,  an  experimental  technique  to  check  the  results  of  the  kine¬ 
matic  rectification  effect  was  devised.  This  technique  involves  applying  a 
circularly  polarized  magnetic  field  transverse  to  the  input  axis  in  order  to 
simulate  the  coning  of  the  gyro  case. 
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IV.  Experimental  Results  Obtained  from  the  Improvement  of  the  Old  Apparatus 

l  .  Modification  of  Apparatus 

Analysis  of  the  data  obtained  from  our  original  apparatus  suggested  that 
the  Tj  measurements  and  possibly  the  free  precession  decay  measurements  were 
strongly  affected  by  a  dipolar  magnetic  field  gradient  arising  from  a  contami¬ 
nant  speck  located  near  the  sample  cell  (see  Figures  5  and  6).  Moreover, 
strong  circumstantial  evidence  indicated  that  this  contaminant  dipole  moment 
was  located  in  the  plastic  pick-up  coil  form.  Since  this  part  of  the  apparatus 
could  be  readily  modified,  it  was  thought  worthwhile  to  replace  all  the  plastic 
components  near  the  sample  cell  with  Macor  parts.  (Macor  is  a  machinable  glass- 
ceramic  which  has  been  found  to  be  magnetically  very  clean.)  The  Macor  field- 
coil  assembly  incorporated  yet  another  improvement  over  the  plastic  assembly 
in  that  it  had  three  orthogonal  field  coils  instead  of  the  previous  two. 

There  were  several  motivations  behind  making  these  improvements  in  the 
original  apparatus.  First,  it  was  hoped  that  the  remaining  apparatus  would 
be  sufficiently  clean  to  allow  it  to  be  operated  as  a  crude  nuclear  gyro. 
Secondly,  if  this  turned  out  to  be  the  case,  to  use  it  to  test  a  number  of 
model  predictions  obtained  in  our  theoretical  analyses.  And  finally,  to  try 
to  use  it  as  a  probe  of  our  new  ultra-low  field  shield.  As  will  be  discussed 
in  the  next  subsection,  it  was  found  that  the  use  of  the  field  coil  assembly 
did  substantially  reduce  the  contaminant-dipole  problem.  This  improvement 
made  it  possible  to  undertake  a  more  extended  and  systematic  study  of  free- 
precession  decay  characteristics  than  had  been  previously  possible.  As  a 
result  of  this  study  it  was  found,  however,  that  the  self-induced  field 
gradient  arising  from  the  interaction  of  the  sample  magnetization  with  the 
third  surface  harmonic  of  the  sample  cell  was  much  larger  than  had  been 
previously  suspected.  The  details  and  implications  of  this  discovery  will 
be  discussed  in  subsection  c. 
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Fig.  6 


Detail  of  sample  cell  and  pickup  coil  form  showing  possible 
location  of  a  ferromagnetic  dipole. 


2.  Measurements  on  Liquid  Mixtures  Obtained  from  the  Modified  Apparatus 

A  sequence  of  longitudinal  relaxation  time  (Tj)  measurements  were  made  on  a 

3  4  3  -4 

standard  He  -He  liquid  mixture  (He  concentration  of  7 x 10  )  as  a  function  of 

the  applied  uniform  field,  Bq,  as  a  diagnostic  check  of  the  modified  apparatus. 

-1  -3/2 

The  results  of  these  measurements  are  plotted  in  Fig.  7  as  Tj  vs.  Bq  .  It 
can  be  shown  that  the  longitudinal  relaxation  rate  of  a  sample  exposed  to  a  di¬ 
polar  perturbation  arising  from  a  macroscopic  moment  located  at  a  distance  from 

the  sample  cell  that  is  small  compared  to  the  cell  diameter  is  of  the  form 

-1  -1  -3/2 

T,  =  T,  +  aB  1  ,  where 
1  1°°  o’ 

a  -  37rm2(yD)  ^  /  (4  /Tvb4). 

In  this  expression,  V  is  the  sample  cell  colume,  D  is  the  He^  diffusion  coeffi- 

3 

cient,  m  is  the  magnitude  of  the  contaminant  dipole  moment,  y  is  the  He  gyro- 
magnetic  ratio,  and  b  is  the  distance  between  the  dipole  and  the  nearest  bound¬ 
ary  of  the  sample  cell.  In  order  for  this  expression  to  be  valid  it  is  necessary 

that  |Sq|  >>  [fj  everywhere  in  the  sample  cell  (B^  is  the  perturbation  field 

2 

arising  from  the  contaminant  dipole),  and  >>  D/yb  . 

If  the  lowest  two  field  points  (taken  at  4.3  and  14  G)  are  omitted  and  the 

data  replotted  on  a  larger  scale  (Fig.  8)  it  is  seen  that  that  data  fits  the 
-3/2 

Bq  dependence  quite  well  for  Bq  42  pG.  A  least  squares  fit  to  this  data 

-12  -1  3/2 

yields  an  estimate  for  the  parameter  a  -3.80  x  10  sec  G  .  This  is  a 

factor  of  five  smaller  than  the  value  obtained  with  the  Delrin  coil  forms  which 

is  also  shown  in  Fig.  8  for  comparison. 

From  this  T^  data  we  can  draw  the  following  conclusions:  1)  The  dominant 

dipolar  perturbation  field  that  was  noted  in  our  original  data  was  indeed  due 

to  magnetic  contaminant  in  the  Delrin  plastic  coil  forms.  2)  The  replacement 

of  the  plastic  coil  forms  has  resulted  in  a  factor  of  five  improvement  in  the 

magnetic-gradient  induced  relaxation  due  to  nearby  dipole  moments.  3)  One  or 

-3/2 

both  of  the  validity  requirements  for  the  Bq  "  dependence  are  violated  in 


Fig.  7 

Fig.  8 


1 


Longitudinal  relaxation  rates  as  a  function  of  the  applied  field 
over  a  range  of  14.2  microgauss  to  1.4  milligauss. 

Same  data  as  in  Fig.  7  plotted  on  a  different  scale.  Lower  curve 
are  the  data  taken  with  Macor  coil  assembly  and  the  upper  curve 
is  a  typical  least-squares  fi*  obtained  with  the  original  plastic 
coil  assembly. 
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fields  less  than  40  pG.  This  last  conclusion  in  conjunction  with  the  measured 
value  of  the  slope  parameter,  a,  leads  us  to  suspect  that  the  residual  dipole 
moment  that  is  still  observed  is  located  either  in  or  on  the  Pyrex  wall  of  the 
sample  cell  itself.  This  suspicion  has  no  direct  implication  for  the  new  pre¬ 
cision  gyro  housing  since  the  new  sample  cell  is  much  larger  and  is  made  of  high 
purity  fused  silica  fabricated  by  precision  grinding  and  lapping  rather  than 
glassblowing.  This  result,  however,  does  point  out  the  importance  of  using  a 
non  magnetic  lapping  compound  (aluminum  oxide  rather  than  rouge  or  cerium 
oxide)  and  rigorous  cleaning  procedures. 

3.  Free  Precession  Decay  Measurement 

a.  Introduction 

Of  more  direct  importance  to  the  gyro  applications  than  the  longitudinal 
relaxation  time  is  the  decay  of  the  transverse  component  of  the  magnetization 
during  free  precession.  It  is  essential  that  long  transverse  decay  times  be 
obtained  in  order  that  precision  rotation-rate  measurements  can  be  made.  In 
what  follows  we  will  first  discuss  the  general  aspects  of  free  precision  decay 
analysis  and  the  kind  of  information  we  hoped  to  obtain  from  this  type  of 
measurement.  Experimental  results  obtained  from  our  modified  apparatus  will 
then  be  discussed  in  the  light  of  two  different  kinds  of  analysis.  These  new 
results  represent  a  significant  improvement  over  our  previous  measurements  in 
that  they  include  some  transverse  decay  times  that  are  an  order  of  magnitude 
longer  than  any  we  had  previously  obtained. 

As  is  discussed  in  detail  in  the  thesis  of  M.  A.  Taber,  there  are  two 

limiting  physical  circumstances  where  free  precession  decay  (fpd)  data  are 

readily  analyzed.  These  two  cases  are  separated  by  the  motional  narrowing 
2  2  2  + 

boundary  y  t  +  1 ,  where  B^  is  the  perturbation  field  responsible  for  trans¬ 
verse  decay  and  tc  is  the  correlation  time  associated  with  the  perturbation. 
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The  perturbation  field  is  a  random  variable  in  time  because  of  the  diffusive 
motion  of  a  given  nuclear  spin  throughout  the  sample  cell  and  is  defined  to  have 

■f 

a  zero  mean.  Bj  can  be  either  microscopic  or  macroscopic  in  origin.  In  the 
former  case,  the  perturbation  field  arises  from  a  nearby  nuclear  or  electronic 
spin,  and  in  a  rough  sense  X^  is  on  the  order  of  the  time  necessary  to  diffuse 
an  atomic  diameter.  In  the  latter  case,  is  the  inhomogeneous  component  of 
the  macroscopic  field  and  x^  is  roughly  the  time  necessary  for  an  atom  to  diffuse 
from  one  side  of  the  sample  cell  to  the  other.  In  our  situation,  the  effect  of 
microscopic  interactions  on  the  transverse  decay  is  quite  weak  compared  to  that 
of  macroscopic  field  gradient,  and  the  current  discussion  deals  exclusively  with 
the  macroscopic  case. 

2  2  2 

When  the  motional  narrowing  condition  y  Bj  xc  <<  1  is  well  satisfied, 
motion  of  the  individual  He ^  atoms  is  sufficiently  rapid  with  respect  to  the 
rate  that  the  perturbation  field  can  reorient  a  nuclear  spin  that  the  perturba¬ 
tion  field  is  largely  averaged  out  and  has  only  a  second  order  effect.  In  this 
limit  the  sample  magnetization  is  uniform  throughout  the  sample  and  the  trans¬ 
verse  decay  can  be  analyzed  by  standard  density  matrix  (Redfield)  spin  relaxa¬ 
tion  theory.  The  theory  yields  a  time  dependence  of  the  transverse  component 
of  the  magnetization  of  the  form  exp(-t/T,r) .  The  mot i onal lv -narrowed  limit 

*-(j 

represents  the  most  desirable  circumstance  with  regard  to  nuclear  gyro  applica¬ 
tions  in  that  the  only  effect  that  the  perturbation  (gradient)  field  has  on 

the  1, armor  frequency  is  a  second  order  frequency  shift  that  is  constant  in  time 

4 

if  the  gradient  itself  is  time- independent . 

If  the  motional  narrowing  condition  is  strongly  violated,  the  fpd  has  a 
significantly  different  character.  In  this  case  the  diffusion  is  relatively 
slow  and  the  spins  process  with  a  l.armor  frequency  given  by  the  local  field 
magnitude.  The  sample  magnetization  then  becomes  nonuniform.  In  this  limit 
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the  effect  of  spin  diffusion  on  the  fpd  is  secondary.  If  it  is  assumed  that  all 
portions  of  the  sample  cell  are  equally  coupled  to  the  pickup  coil  (which  is  a 
good  approximation  with  the  modified  coil  assembly)  and  the  gradient  field  is 
such  that  ! Bj (r) j  <<  Bq  everywhere  in  the  sample  cell,  then  the  transverse  fpd 
signal  represents  the  Fourier  transform  of  the  lineshape  function.  The  line- 
shape  function  is  an  unnormalized  distribution  function  of  the  z  component  of 

-v  -*■  ->  •+ 

B  +  B.  (r)  over  the  sample  volume  where  B  is  taken  to  be  in  the  z  direction. 

An  fpd  arising  from  this  kind  of  situation  is  not  as  useful  for  nuclear  gyro 
applications  as  a  motional ly-narrowed  fpd,  since  any  lineshape  asymmetry  can 
cause  a  time-dependent  phase  shift  which  would  simulate  a  rotation.  In  spite 
of  this,  however,  an  fpd  obtained  from  a  situation  where  the  motional  narrowing 
criterion  is  strongly  violated  does  not  represent  useless  information,  since  it 
can  be  readily  analyzed  to  provide  useful  diagnostic  information;  i.e..  the 
field  distribution  function. 


Because  of  the  importance  of  the  motional  narrowing  criterion,  it  is 

useful  to  make  things  more  concrete  by  specifying  some  of  the  numbers  which 

typify  our  experimental  situation.  In  the  case  of  a  spherical  geometry,  the 

2 

slowest  diffusion  mode  has  a  time  constant  r..  =  0.2308  R  /D,  where  R  is  the 

11  o  ’  o 

radius  of  the  sample  cell  and  D  is  the  He  diffusion  coefficient.  In  our  case 

-4  2-1  3  4 

Rq  =  0.5  cm  and  D  =  10  cm  sec  for  He  in  liquid  He  .  Thus  the  longest 
correlation  time  associated  with  diffusive  motion  is  x  ^  =  10  min.  If  we  now 
consider  for  the  sake  of  definiteness  a  uniform  gradient  of  the  form  (r)  = 

^  A  a  2-*.  2  2 

g(ix  -  Sjjy  -  ^kZ),  then  the  motional  narrowing  criterion  becomes  y  Tjl  = 
222  2  2262 

(3/10)Ro  g  y  xn =  0.02  g  y  R  D  «1.  For  the  numbers  that  have  been  speci- 

-7  -1 

fied  this  means  that  g«  3  x  10  G  cm  .  If  it  is  assumed  that  the  motional 


narrowing  condition  is  indeed  satisfied,  then  an  expression  for  the  transverse 
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decay  time  constant  can  also  be  obtained  for  our  hypothetical  gradient: 
T2g_1  =  JsC0.0229)Y2g2RQ4D'1  , 

where  it  has  been  very  reasonably  assumed  that  >>  1 .  If  this  expression 

is  combined  with  the  motional  narrowing  criterion,  then  it  is  seen  that  T  >> 

2G 

6  Tj^.  Thus  the  signature  of  motional  narrowing  in  our  experimental  environ¬ 
ment  would  be  an  fpd  with  an  exponential  envelope  having  a  time  constant  of 
more  than  one  hour. 

b.  Fourier  Analysis  of  Experimental  FPD  Data  on  a  Liquid  Mixture 

3  4  3 

Free  precession  data  of  a  standard  He  -He  liquid  mixture  (He  concentration 
-4 

ot  7  x  10  )  were  initially  taken  on  a  chart  recorder  but  were  subsequently  taken 

by  analog-to-digital  conversion  in  order  to  facilitate  discrete  Fourier  analysis 
by  computer.  The  technique  used  was  to  start  with  the  magnetization  initially 
aligned  with  the  applied  field,  The  applied  field  was  then  adiabatical ly 

shifted  in  direction  by  an  angle  equal  to  the  desired  precession  cone  half-angle 
followed  by  suddenly  restoring  B^  to  its  original  direction.  The  magnetization 
then  precessed  in  a  cone  about  Bq  until  the  transverse  component  decayed  to  zero. 
The  precession  core  half-angle  was  selected  to  be  45°  since  this  gave  a  signal 
that  was  70%  of  the  maximum  attainable  (with  a  half-angle  of  90°)  but  resulted 
in  the  loss  of  only  30%  of  the  magnetization  per  each  fpd  measurement .  The 
loss  in  the  longitudinal  component  of  the  magnetization  was  small  since  the 
longitudinal  relaxation  time  (T^ )  was  typically  much  longer  than  even  the  longest 
fpd  measurements.  In  all  the  measurements  to  be  discussed  here,  the  value  of 

i  -V  , 

B  |  during  the  fpd  was  ^elected  to  be  61.7  pG  so  that  the  Larmor  frequency  was 
0.2  Hz. 


36 


The  digitally  recorded  data  were  analyzed  using  standard  discrete  Fourier 
transform  (dft)  techniques.  Because  of  the  one-sided  nature  of  the  data 
record  (i.e.,  no  data  exists  before  t  =  0) ,  however,  some  care  had  to  be 
taken  to  remove  phase  shifts  and  time  delays  before  symmetrizing  the  data 
and  performing  the  dft.  Data  which  were  obtained  from  a  sequence  of  fpd 
measurements  made  during  a  single  run  and  were  thus  analyzed  are  shown  in 
Figs.  9-15.  There  are  only  two  parameters  which  vary  between  these  fpd  trans¬ 
forms:  the  direction  of  the  applied  uniform  field,  B  ,  and  the  magnitude  of 
the  magnetization.  In  the  first  and  second  fpd,  Bq  was  in  the  vertical  direc¬ 
tion  (i.e.,  the  fill  line  direction),  and  in  all  other  fpd's,  Bq  was  in  the 
horizontal  direction.  The  magnitude  of  the  magnetization  also  monotoni cal  2 y 
declines  from  one  fpd  data  record  to  the  next  because  of  the  loss  of  the  trans¬ 
verse  component. 
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Figs.  9-15.  Fourier  transforms  of  free  precession  decays  (lineshape  functions) 
of  a  liquid  mixture  as  a  function  of  initial  BmQ  =  8ttM/3,  where  M 
is  the  average  sample  magnetization.  Magnetometer  sensing  direc¬ 
tion  is  horizontal  and  is  designated  the  x  direction.  The  z  direc¬ 
tion  is  taken  to  be  vertically  upward.  In  all  cases,  the  static 


field  magnitude  is  Bq  =  01.7  pG. 
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Fig.  12 
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Fig.  13 
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Fig.  14 


There  are  a  number  of  noteworthy  aspects  of  these  1 inshape  curves.  First, 
the  linewidths  decreased  markedly  as  the  magnitude  of  the  magnetization  decreased. 
This  effect  was  noted  during  the  fpd  data  taking  and  forced  longer  and  longer 
(up  to  40  minutes)  data  records  to  be  made.  Secondly,  there  are  features  of 

the  lineshape  curves  which  are  physically  implausible  if  they  are  viewed  in  a 

*+  ->  -+■ 

straightforward  manner  as  distribution  functions  of  the  z  component  of  Bq  +  B^ (r) . 
In  particular,  the  curves  have  "wings"  where  the  lineshape  function  changes  signs, 
and  in  addition  there  frequently  is  a  very  sharp  jump  in  the  middle  of  the  line. 

It  was  at  first  thought  that  these  features  were  artifacts  due  to  errors  in  the 
analysis.  Subsequent  work,  however,  made  it  clear  that  the  analysis  was  correct 
but  that  the  simple  distribution  function  model  was  totally  inadequate. 

The  reason  for  the  inadequacy  of  the  simple  model  is  that  the  large  scale 
gradient  (as  opposed  to  the  localized  dipolar  perturbation  which  dominates  the 
Tj  data)  is  dominated  by  the  gradient  arising  from  the  sample  magnetization  due 
to  asphericity  of  the  sample  cell  shape.  In  what  follows  we  will  describe  in 
somewhat  more  detail  the  mechanism  for  the  generation  of  the  gradient.  We  will 
then  discuss  in  qualitative  terms  how  a  gradient  arising  from  the  magnetization 
itself  introduces  significant  complexities  beyond  the  normal  inhomogeneous  line 
broadening  that  is  usually  encountered  in  spin  resonance  experiments.  The 
problem  is  sufficiently  complex  that  a  more  quantitative  analysis  has  not  been 
undertaken  since  the  new  precision  sample  cell  will  reduce  the  effect  by  many 
orders  of  magnitude. 

c.  Self-Induced  Broadening  Due  to  Sample-Shape  Asphericity 

The  internal  macroscopic  B-field  produced  by  a  uniformly  magnetized  sample 
is  in  general  not  uniform  unless  the  sample  shape  is  ellipsoid,  spheroid  or  a 
sphere.  In  most  spin  resonance  experiments  this  is  of  little  consequence  since 
the  sample  magnetization  is  sufficiently  small  that  the  gradient  thus  produced 
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is  trivial  compared  to  the  gradient  in  the  applied  field.  In  our  case,  however 
the  combination  of  a  large  magnetization,  very  small  gradients  in  the  applied 
field,  and  the  sample  shape  have  conspired  to  make  the  opposite  true. 


In  the  thesis  of  M.  A.  Taber,  a  calculation  is  outlined  which  allows  the 
first-order  estimation  of  the  internal  magnetic-field  structure  of  a  nearly 
spherical,  uniformly  magnetized  sample.  In  general  this  calculation  shows  that 
if  the  deviation  of  the  sample  shape  from  spherical  is  expanded  in  spherical 
harmonics  and  if  the  internal  magnetic  field  is  specified  in  terms  of  a  magneti 
scalar  potential  which  is  likewise  expanded  in  spherical  harmonics,  then  the 
coefficients  of  the  two  expansions  can  be  related  to  each  other  in  first  order 
by  the  usual  quantum  mechanical  angular  momentum  coupling  coefficients.  The 
third  "angular  momentum"  turns  out  to  have  the  value  £  =  1.  Thus,  by  using 
triangle  rule  and  parity  consideration  one  can  quickly  obtain  an  intuitive 
assessment  of  which  sample  cell  shapes  will  be  important  in  generating  a  speci¬ 
fied  field  structure.  In  our  particular  case,  we  are  mostly  interested  in  the 
uniform  gradient  (Y2m^  terms  in  the  field  structure.  If  the  coordinate  system 
is  chosen  such  that  the  coefficients  of  the  terms  in  the  cell-shape  expan¬ 
sion  are  zero,  the  only  terms  that  have  a  first-order  effect  in  producing  a 
uniform  gradient  are  the  Y^m  terms.  Since  Y^q  represents  a  pear-shape  pertur¬ 
bation,  it  is  clear  from  the  shape  of  our  sample  cell  (Fig.  6)  that  the  Y^m 
terms  are  likely  to  be  very  important  if  not  totally  dominant  in  the  cell-shape 
expansion.  Thus,  it  is  not  too  surprising  that  a  gradient  that  is  proportional 
to  the  sample  magnetization  has  been  observed  in  our  apparatus. 

The  fact  that  the  magnetic  field  gradient  is  related  to  the  sample  magneti 
zation  has  a  number  of  significant  implications.  The  most  obvious  one  is  that 
as  the  magnitude  of  the  sample  magnetization  decreases,  the  linewidths  would 


get  narrower  as  we  have  experimentally  observed.  By  writing  the  magnetization 
as  M  =  Mj  +  fS^,  where  Mj  is  the  component  of  M  in  the  direction  of  ?  and 
is  the  precessing  transverse  component,  a  number  of  other  implications  can  be 
perceived.  First,  as  decays  during  a  fpd  the  the  magnitude  of  the  gradient 
can  be  expected  to  decrease  and  its  structure  change.  Qualitatively  one  would 
expect  the  fpd  to  show  a  rather  rapid  initial  dephasing  followed  by  a  long  "tail" 
during  which  the  decay  is  much  slower.  This  type  of  behavior  can  explain  the 
presence  of  the  sharp  dispersive  feature,  which  arises  from  the  long  tail,  in 
the  middle  of  a  broader  line,  which  is  due  to  the  initial  rapid  decay.  (The 
fact  that  the  sharp  feature  is  a  jump  rather  than  a  peak  is  due  to  the  fact  that 
the  phase  of  the  tail  is  shifted  with  respect  to  the  initial  precession  phase.) 
The  presence  of  the  negative  wings,  however,  remains  obscure. 

Another  implication  of  the  self-induced  broadening  model  is  with  regard 
to  the  spatial  uniformity  of  the  transverse  component  of  the  magnetization. 
Initially  the  magnetization  would  start  out  fairly  uniform,  but  as  the  preces- 

4 

sion  proceeds,  dephasing  will  cause  significant  inhomogeneities  in  to  arise. 

At  this  point,  of  course,  it  is  no  longer  possible  to  simply  estimate  the  inter¬ 
nal  macroscopic  field  structure.  The  presence  of  spin  diffusion,  however,  tends 

-4 

to  reduce  the  spatial  variation  of  M^,  and  at  the  same  time  causes  the  dephasing 
to  be  irreversible.  If  is  not  too  large  (or  is  optimally  oriented  with 

-4 

respect  to  the  sample  cell  shape)  and  the  average  of  becomes  sufficiently 
small,  it  may  become  possible  for  the  internal  field  gradient  to  be  small  enough 
to  approach  the  motional  narrowing  condition.  In  this  case,  the  magnetization 
would  once  again  become  fairly  uniform  and  the  tail  of  the  fpd  would  be  expo¬ 
nential  with  a  long  delay  time  with  respect  to  the  diffusion  time.  Because  our 
data  in  some  of  the  latter  field  measurements  had  an  appearance  that  suggested 
that  the  motional  narrowing  condition  was  at  least  being  approached,  a  separate 


] 
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type  of  analysis  was  carried  out  to  investigate  this  possibility. 

d.  Envelope  Analysis  of  FPD  Tails 

The  raw  data  from  the  last  four  fpd's  where  the  initial  magnitude  of  the 
internal  B-field  due  to  the  sample  magnetization  ranged  between  14  and  4  pG  was 
subjected  to  a  digital  detection  and  filtration  routine  which  yielded  the  ampli¬ 
tude  (envelope)  of  the  fpd  signal  as  a  function  of  time.  The  logarithm  cf  the 
amplitude  was  then  plotted  and  inspected  for  linear  behavior  in  the  latter  por¬ 
tion  of  the  record.  In  three  of  the  four  records  which  were  examined  this  way 
it  was  found  that  the  amplitude  of  fpd  tail  did  indeed  decay  exponentially.  In 
one  exceptional  case  the  onset  of  the  exponential  tail  was  anomalously  delayed, 
and  since  the  sequence  of  events  prior  to  recording  this  particular  fpd  was 
significantly  different  from  the  others,  this  particular  record  was  omitted 
from  further  analysis.  The  remaining  three,  however,  were  subjected  to  a  least 
squares  fit  to  an  exponential  decay.  An  example  of  one  of  the  fit.,  is  shown  in 
Fig.  16  and  in  Fig.  17  the  value  of  the  three  transverse  decay  rates  (T2  *) 
thus  obtained  are  plotted  as  a  function  of  the  magnitude  of  the  initial  average 
internal  B-field  due  to  the  sample  magnetization  (B^  =  8ttM/3). 

Fig.  17  indicates  that  the  transverse  decay  rate  is  a  monotonic  function 
of  B .  This  is  presumably  due  to  the  self-induced  gradient  still  being  a  signi¬ 
ficant  contributor  to  the  relaxation  rate.  The  longest  T^  was  21  minutes,  which 
was  obtained  for  B  =4  yG.  Since  the  slowest  diffusion  mode  decay  time  is 
approximately  10  minutes,  it  appears  that  this  represents  being  on  the  border 
of  motional  narrowing.  It  will  be  recalled  that  it  would  be  necessary  for  T0 
to  be  in  excess  of  an  hour  in  order  for  motional  narrowing  to  be  clearly  esta¬ 
blished.  Another  indication  of  the  marginal  nature  of  the  situation  is  the  fact 
that  one  data  record  (the  rejected  one)  was  adversely  affected  by  events  occur¬ 
ring  before  the  fpd  was  initiated.  This  is  a  good  example  of  how  difficult  it 
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Fig.  17.  Plot  of  the  natural  logarithm  of  the  envelope  of  a  freee  precession 
decay  tail  together  with  the  least-squares  fit.  In  this  particular 
case  initial  Bmo  =  3.97  microgauss,  and  =  23  min. 

Fig.  18.  T?  of  the  free  precession  decay  tails  for  three  different  values  of 
the  initial  B 

mo 
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is  to  obtain  understandable  transverse  decay  results  from  liquid  mixtures  where 

inhomogeneities  in  the  sample  magnetization  take  a  long  time  to  dissipate. 

e.  Measurements  made  on  a  Gas  Mixture  Sample 

3  4 

There  were  two  reasons  to  try  using  a  He  -He  gas  mixture  despite  the 

known  deficiencies  of  the  sample  cell.  The  first  is  that  the  He3  diffusion 

coefficient  could  be  expected  to  be  an  order  of  magnitude  or  more  larger  in  the 

4 

gas  than  in  the  liquid.  (The  ratio  of  He  liquid  density  at  4.2  K  to  that  of 
the  saturated  vapor  is  approximately  7.5).  Thus  the  use  of  a  gas  mixture  should 
help  answer  the  question  of  whether  the  liquid  mixtures  were  approaching  motional 

4 

narrowing.  The  second  reason  was  to  determine  whether  the  reduced  amount  of  He 
would  increase  the  loss  of  He3  polarization  during  the  sample  cell  filling  pro¬ 
cess.  This  issue  is  of  interest  because  the  new  precision  gyro  apparatus  is 
designed  to  be  used  with  gas  mixtures.  Because  of  this  latter  issue,  the  amount 

4 

of  He  chosen  for  these  measurements  was  that  necessary  to  produce  a  final  pres¬ 
sure  slightly  less  than  one  atmosphere.  This  is  the  maximum  amount  of  He4  that 
can  be  used  without  forming  liquid  in  the  sample  cell.  It  should  also  be  noted 

that  the  He3  density  will  be  lower  in  a  gas  mixture  than  in  a  liquid  mixture. 

3  4 

The  ratio  of  He  concentrat ion  in  He  vapor  at  4  K  to  that  in  the  liquid  is  1.4 
to  1.  This  effect,  however,  is  overwhelmed  by  the  much  greater  density  of  the 
liquid  compared  to  that  of  the  vapor  with  the  net  result  that  the  He3  density 
in  the  liquid  is  approximately  five  times  that  in  the  gas. 

We  now  summarize  the  results  that  were  anticipated  in  the  gas  mixture 
measurements : 

1)  The  magnetization  should  be  reduced  to  30%  of  the  liquid  value  (every¬ 
thing  else  being  equal)  because  of  the  lower  He3  density  if  the  reduced  amount 

4  3 

of  He  does  not  adversely  affect  the  relaxation  of  He  during  the  sample  cell 

filling  process. 
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2)  might  well  decrease  because  the  reduced  He  density  could  increase 
the  efficiency  of  wall-induced  relaxation.  There  is  insufficient  information  to 
judge  how  the  gradient-induced  relaxation  rate  would  change. 

3)  Despite  a  possible  decrease  in  it  is  expected  that  the  transverse 
decay  time  should  increase  for  two  reasons:  a)  reduction  in  magnetization  will 
yield  a  smaller  cell-shape- induced  gradient,  and  b)  the  more  rapid  diffusion  will 
cause  motional  narrowing. 

In  one  experimental  run  we  were  able  to  make  one  Tj  measurement  at  61.7  uG 
and  a  sequence  of  free  precession  decay  measurements  in  exactly  the  same  fashion 
described  for  the  liquid  mixture.  The  initial  magnetization  was  reasonably 
within  expectation.  The  initial  magnetization  yielded  -  9  pG  which  is  about 
20%  of  the  value  obtained  in  a  fairly  comparable  liquid  mixture  run.  Since  we 

4 

had  expected  a  30%  factor,  this  is  an  indication  that  the  reduced  amount  of  He 
may  have  led  to  a  somewhat  less  efficient  transfer  of  polarized  He^  from  the 
optical  pumping  cell  to  the  sample  cell.  This  effect,  however,  does  not  appear 
to  be  large  enough  to  cause  concern  about  the  viability  of  future  gas  mixture 
work  using  .hese  techniques. 

The  Tj  measurement  yielded  a  relaxation  time  of  2.7  hours.  This  is  a  factor 
of  seven  less  than  would  be  expected  for  a  liquid  sample  in  a  field  of  61.7  ,.G. 

As  was  the  case  with  the  liquid  mixture,  it  is  fairly  certain  that  the  Tj  was 
dominated  by  surface  effects,  that  is,  wall -induced  relaxation  and  superficial 
magnetic  dipole  gradients.  Thus  as  we  have  previously  noted,  the  T^  in  the  pre¬ 
cision  gyro  apparatus  should  be  significantly  longer  because  a  much  larger  cell 
with  a  solid  H  coating  will  be  used. 

Although  the  magnetization  and  Tj  measurements  in  the  gas  mixture  were  not 
very  surprising,  this  was  not  true  of  the  transverse  decay  measurements.  These 
measurements  were  made  and  analyzed  using  virtually  the  same  techniques  as  were 
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used  before.  Both  discrete  Fourier  transform  and  envelope  decay  analysis  are 
shown  in  Figs.  18-23.  The  most  remarkable  of  these  lineshape  plots  are  those 


obtained  with  Bq  in  the  z  direction  (Figs.  18,  19).  These  two  lineshapes  are 
distinctly  bimodal  in  character,  suggesting  that  there  were  two  separate  regions 
in  the  cell  with  different  Larmor  frequencies.  Note,  however,  that  the  remain¬ 
ing  lineshapes  which  resulted  when  Bq  was  in  a  horizontal  direction  are  very 
narrow  with  little  or  no  indication  of  a  secondary  or  satellite  peak.  The  only 
exception  to  this  was  the  first  measurement  made  with  a  horizontal  field  (Fig. 

20)  which  shows  a  slight  shoulder  indicating  the  possibility  of  an  unresolved 
secondary  peak. 

There  appears  to  be  only  one  straightforward  explanation  for  these  somewhat 

4 

peculiar  results.  That  is  that  there  was  slightly  too  much  He  ,  and  a  small 
fraction  of  the  sample  cell  (the  data  suggest  ~ 10%)  contained  liquid  mixture 
rather  than  gas.  The  liquid  portion  would  have  a  magnetization  approximately 
five  times  that  of  the  gas  part  and  would  have  the  shape  of  a  rather  aspheric 
puddle  in  the  bottom  of  the  sample  cell.  These  two  factors  mean  that  the  Larmor 
frequency  of  the  liquid  portion  of  the  sample  cell  would  be  shifted  from  that  of 
the  gas  part  by  an  amount  proportional  to  the  average  sample  magnetization.  (See 
ref.  4  for  a  detailed  discussion  of  the  effect  of  sample  shape  on  the  Larmor 
frequency).  This  is  exactly  what  was  observed  in  the  vertical  field  measure¬ 
ments  where  a  34%  decrease  in  the  average  sample  magnetization  was  accompanied 
by  a  32%  decrease  in  the  splitting  between  the  peaks.  (The  fact  that  the  satel¬ 
lite  peak  shifted  from  one  side  of  the  primary  peak  to  the  other  is  consistent 
with  the  fact  that  Bq  was  reversed  between  measurements  while  the  magnetization 
was  not.) 
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Figs.  18-23.  Fourier  transforms  of  free  precession  decays  (lineshape  functions) 
of  a  nominal  gas  mixture  (see  text)  as  a  function  of  initial 
B  ^  =  8irM/3  where  M  is  the  average  sample  magnetization. 
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An  attempt  to  model  this  proposed  explanation  by  considering  the  liquid 

portion  of  the  sample  to  be  an  oblate  spheroid  located  adjacent  to  a  spherical 

gas  sample  yielded  qualitatively  correct  results.  Quantitatively,  however,  the 

agreement  was  not  totally  convincing.  In  particular  the  model  predicted  a  '-40% 

— ► 

larger  splitting  than  that  observed  for  a  vertical  Bq.  Additionally,  although 
the  model  correctly  indicates  that  the  splitting  should  be  smaller  for  a  hori¬ 
zontal  B  ,  it  also  indicates  that  it  should  be  readily  resolvable  for  the  larger 
magnetization  values.  It  should  be  noted,  however,  that  this  model  is  rather 
crude  and  has  at  least  one  major  defect:  the  exchange  of  He"5  atoms  between  the 
two  sample  regions  has  been  entirely  ignored.  The  effect  of  such  exchange 
would  be  to  tend  to  fuse  the  two  peaks  into  one  broader  one. 

Thus,  while  the  identification  of  the  bimodel  lineshape  of  the  vertical 
precession  data  with  a  liquid-gas  phase  separation  is  not  definitive,  there 
is  sufficient  reason  to  be  cautious  in  drawing  conclusions  from  this  data.  In 
spite  of  this  limitation,  an  envelope  decay  analysis  was  performed  on  all  of  the 
free  precession  decay  measurements  made  with  a  horizontal  B^.  Of  these  the  last 
three  (corresponding  to  Figs.  21-23)  were  accurately  fit  by  a  simple  exponential 
decay.  If  in  fact  there  were  liquid-gas  phase  separation,  this  would  indicate 
that  for  these  weak  magnetizations  the  exchange  of  He^  atoms  between  liquid  and 
vapor  was  sufficiently  rapid  to  completely  wash  out  the  difference  in  Larmor 
frequencies  between  the  two  phases.  In  any  case,  it  was  found  that  the  exponen¬ 
tial  decay  rate  (T^  varied  as  the  square  of  the  mean  sample  magnetization 
with  the  longest  measured  time  constant  being  33  minutes.  This  is  exactly  the 
kind  of  dependence  that  would  be  expected  if  the  transverse  decay  were  controlled 
by  the  rapid  exchange  of  spins  between  two  regions  of  slightly  differing  Larmor 
frequencies.  It  is  also,  however,  what  would  be  expected  if  the  transverse 
relaxation  were  dominated  by  a  field  gradient  arising  from  sample-shape  aspheri- 
city . 
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Because  of  the  possibility  of  liquid-gas  phase  separation  it  is  not  poss¬ 
ible  to  draw  any  firm  conclusions  by  making  comparisons  with  the  liquid-mixture 
transverse  decay  times.  It  should  be  noted,  nevertheless,  that  with  the  excep¬ 
tion  of  the  33  min.  decay  time  (the  longest  that  we  have  observed)  the  transverse 
decay  times  of  the  gas  (or  gas/liquid''  mixture  were  shorter  than  those  of  the 
liquid  mixture,  despite  the  fact  that  all  the  magnetizations  used  in  the  former 
measurements  were  less  than  those  in  the  latter. 

V.  Summary  and  Conclusions 

In  summary,  progress  has  been  achieved  in  the  following  four  areas: 

a)  A  new  ultra-low  magnetic  field  shield  has  been  established  which 

__  g 

achieved  a  field  of  2  x  10  G  over  a  substantial  volume.  This  not  only  exceeds 
our  requirements,  but  to  the  best  of  our  knowledge  is  the  lowest  absolute  magne¬ 
tic  field  region  ever  made. 

b)  The  precision  fused  quartz  gyro  housing  has  been  received  and  is  ready 
for  final  tumble-lapping,  shield  deposition,  and  assembly  procedures.  The  new 
cryostat  probe  is  in  the  construction  phase. 

c)  The  engineering  Ph.D.  thesis  of  Capt.  Gerald  Shaw  which  analyzes  several 
aspects  of  the  dynamic  behavior  of  our  nuclear  gyro  has  been  completed. 

d)  Our  original  experimental  apparatus  was  modified  and  run  four  more 
times  (once  with  a  nominal  gas  mixture)  with  the  result  of  achieving  an  exponen¬ 
tial  decay  characteristic  in  the  limit  of  a  small  magnetization  and  a  signifi¬ 
cant  reduction  in  the  magnetic  contaminant  problem.  This  has  had  the  effect  of 
increasing  the  transverse  decay  time  by  over  an  order  of  magnitude  over  what  we 
had  previously  observed.  The  data  indicated  that  the  transverse  decay  times 
were  strongly  affected,  if  not  dominated  by  sample-induced  gradients. 


The  improvement  achieved  with  the  modification  of  our  original  apparatus 
has  made  clear  the  impact  of  the  sample  shape  on  the  transverse  decay  time  and 
has  convinced  us  of  the  importance  of  using  gas  mixtures  in  a  precision  sample 
cell  made  out  of  high  grade  material.  Since  these  are  features  of  the  new 
apparatus,  we  feel  that  if  the  blown  Pyrex  cell  worked  as  well  as  it  did,  we 
have  every  expectation  that  the  new  apparatus  will  achieve  the  conditions 
necessary  for  quality  gyro  performance. 
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